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To understand the role that degeneracy, hybridization, and nesting play in the magnetic and 
pairing properties of multiorbital Hubbard models we here study numerically two types of two- 
orbital models, both with hole- like and electron-like Fermi surfaces (FS's) that are related by nesting 
vectors (tt, 0) and (0, n). In one case the bands that determine the FS's arise from strongly hybridized 
degenerate d^j^ and dy^ orbitals, while in the other the two bands are determined by non-degenerate 
and non-hybridized s-like orbitals. Using a variety of techniques, in the weak coupling regime 
it is shown that only the model with hybridized bands develops metallic magnetic order, while 
the other model exhibits an ordered excitonic orbital-transverse spin state that is insulating and 
does not have a local magnetization. However, both models display similar insulating magnetic 
stripe ordering in the strong coupling limit. These results indicate that nesting is a necessary 
but not sufficient condition for the development of ordered states with finite local magnetization in 
multiorbital Hubbard systems; the additional ingredient appears to be that the nested portions of the 
bands need to have the same orbital flavor. This condition can be achieved via strong hybridization 
of the orbitals in weak coupling or via the FS reconstruction induced by the Coulomb interactions 
in the strong coupling regime. This effect also impacts the pairing symmetry as demonstrated by 
the study of the dominant pairing channels for the two models. 
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I. INTRODUCTION 



Among the several aspects of the study of the iron- 
based superconductors that are still controversial and un- 
settledji the following two questions have attracted con- 
siderable attention: (i) Does the magnetic order observed 
in the parent compounds^ arise from the nesting proper- 
ties of the non-interacting (or high temperature) Fermi 
surface^i^ or should a better description be based on the 
superexchange Heisenberg interactions between localized 
magnetic moments?— (ii) What is the pairing mechanism, 
to what extent is the pairing symmetry determined by 
nesting and, what is the actual symmetry and momen- 
tum dependence of the pairing operator? In particular, 
what is the role that the orbital degrees of freedom play 
in this context? 

The origin of the magnetic state is being vigorously 
debated. One proposal, based on fermiology, is the exci- 
tonic mechanism in which electron-hole pairs arc formed 
by one electron and one hole from different FS's nested 
with nesting vector Q. In this context some studies disre- 
gard the orbital structure of the bands^i^"— while others 
stress the role played by their orbital compositioni^"— 
Another approach focuses on the order of the localized 
moments that develop in the presence of strong Coulomb 
interaction o^'^^" — and relies on ab initio result a^^i^° that 
suggest that the pnictides are moderately, rather than 
weakly, correlated, conclusion supported by photoemis- 



sion measurements indicating mass enhancements due to 
electron correlations as large as 2-3 »^ 

The pairing mechanism in the pnictides is also con- 
troversial. Most of the pairing operators that have been 
proposed in the literature either ignore the multi-orbital 
characteristics of the problem or consider Cooper pairs 
that arc made out of electrons located at the same or- 
bital. A majority of these previous studies have been 
performed in the weak coupling limit. The original pro- 
posal of the s± pairing state dealt with the overall sym- 
metry of the pairing operator but without distinguishing 
among the spatial vs. orbital contributions to its par- 
ticular formi^"^ Other authors^ have considered a spin- 
fiuctuation-induccd pairing interaction and also assumed 
that Cooper pairs are predominantly made of electrons 
in the same orbital. A Random-Phase Approximation 
(RPA) analysisii concluded that the pairing is, again, 
intraorbital, both for the Aig (s-wave) and Big (d-wave) 
symmetries. Among the authors that have used the con- 
ceptually different strong coupling approach, some have 
studied effective single orbital modelsS^ while others in- 
corporated two orbitals^ but still only considering intra- 
orbital pairing operators. The same model was also stud- 
ied under a mean- field approximation^l with the assump- 
tion that exchange takes place between spins on the same 
orbitals and, again, only intraorbital pairs were proposed. 

Among the early first studies of multiband supercon- 
ductors, Suhl et ali^ considered two tight-binding bands. 



hypothctically identified with s and d orbitals, and the 
effect of weak electron-phonon interactions. Under these 
assumptions, it was reasonable to expect that the Cooper 
pairs would be formed by electrons belonging to the same 
band. However, the actual orbital composition of the 
pairs was not addressed. The interacting portion of the 
Hamiltonian was written in the band representation and 
this model was proposed by analogy with models used 
in the BCS theory, assuming that emission and absorp- 
tion of a phonon could occur in four ways. These four 
processes corresponded to pair scattering within each of 
the two FS's and pair hopping from one FS to the other. 
This last process would occur if the exchanged phonon 
has enough momentum to allow the Cooper pair to jump 
from a FS to the other, and it can occur even if the or- 
bitals do not hybridize to form the bandsj^ In this case, 
the expected pairing operator is the traditional on-site 
s-wave state of the BCS theory, with a momentum inde- 
pendent gap. In principle, independent gaps may arise 
on the different FS's2^ unless the orbitals are hybridized 
by the symmetries of the Hamiltonian, in which case the 
gaps will have to be related to each other and obey the 
symmetries of the systemJ^ 

The previous discussion applies to superconductors 
driven by the electron-phonon interaction. However, it is 
believed that the most relevant interactions in the pnic- 
tides are the Coulomb repulsion and Hund magnetic ex- 
change. These interactions are more easily expressed in 
real space and in the orbital representation. In fact, the 
effective form of the Coulomb interaction in the band 
representation is more complicated than the expression 
provided by Suhl et al^ for the electron-phonon interac- 
tion. In particular, it has been showii^ that a pair hop- 
ping term, such as the one introduced by Suhl et al. oc- 
curs only if the orbitals get hybridized to form the bands. 
If the orbitals are not hybridized this type of term is not 
present in the effective interaction Hamiltonian. In addi- 
tion, when the bands are made of hybridized orbitals, as 
it is the case for the iron pnictides?^ the actual orbital 
structure of the pairs needs to be considered since due to 
the Coulomb repulsion on-site pairing is not expected to 
occur, and the overall symmetry properties of the pairing 
operators may be a function of their spatial and orbital 
components jiSiii^ 

To understand the role that the orbitals play in the 
case of electrons with strongly hybridized bands that are 
interacting via the Coulomb repulsion, as believed to oc- 
cur in the case of the pnictides in the context of the mag- 
netic scenario for superconductivity, in this manuscript 
we present and discuss Lanczos numerical, Hartrce mean- 
field, and RPA studies of two different two-orbital mod- 
els, both displaying identical Fermi surfaces. One of them 
is the well-known and widely used two-orbital model for 
the pnictides^ i^"'^^ based on the two strongly hybridized 
degenerate dxz and dy^ orbitals of iron, while the second 
is a two-band "toy- model" (dubbed the s-model) whose 
bands arise from two non-hybridized, non-degenerate, s- 
likc orbitals that is introduced here for the first time. The 



latter model has a FS qualitatively similar to that of the 
pnictides. In both cases a hole (electron) FS is located at 
the r/M (X/Y) points of the Brillouin zone (BZ). The 
hole and electron FS's are connected by nesting vectors 
(tt, 0) and (0, tt). The role that the nesting and the or- 
bitals play in the magnetic and pairing properties of these 
models will be here investigated and discussed, both in 
the weak and strong coupling regimes. 

Besides its conceptual relevance, the results pre- 
sented here should also be framed in the context of 
recent bulk-sensitive laser angle-resolved photoemission 
(ARPES) experiments^^ on BaFe2(Aso.65Po.M]2 and 
Bao.6Ko.4Fe2As2. The main conclusion of Ref. [2a is the 
existence of orbital independent superconducting gaps 
that are not expected from spin fluctuations and nest- 
ing mechanisms, but are claimed to be better explained 
by magnetism-induced interorbital pairing and/or orbital 
fluctuations. This is argued based on the observation 
that the 3z^r^ orbital that forms one of the hole pock- 
ets at the BZ center, but that docs not have a nested 
partner with the same orbital at the electron pock- 
ets, nevertheless appears to develop a superconducting 
gap. Another interesting experimental result that chal- 
lenges the role of nesting in the physics of the pnic- 
tides is a careful measurement of the de Haas- van Alplicn 
(dHvA) effect in BaFe2P2, the end member of the series 
BaFe2(Asi_a;Pa:)2, indicating that this non-magnetic and 
non-superconducting compound displays the best nesting 
of all the compounds in the scries i^ 

The manuscript is organized as follows. In Section |ll] 
the models are introduced. The magnetic properties are 
presented in Section IIIII while the pairing properties are 
the subject of Section Hvl Section |V] is devoted to the 
conclusions. 



II. MODELS 



A. d-model 

The reference model that will be considered here is the 
widely-used two-orbital model^ '^^i^^ based on the d^z {x) 
and dyz [y) Fe orbitals of the pnictides. Since the two 
orbitals are degenerate, an important detail is that the 
direction along which each orbital is defined is actually 
arbitrary. Two directions have been used in the litera- 
ture: X, ySiiiSii^ with the x and y axes along the directions 
that connect nearest-neighbor iron atoms, and AT, F— 12£ 
with the X and Y axis rotated 45° with respect to the 
(x,2/) set. In terms of the d^^r and dyz orbitals the tight- 
binding dispersion of the two-orbital model is given by^ 



£,xyi}i) = [-{h + i2)(cos kx + cos fc^)- 

4^3 cos kx cos ky — ^]tq 

— (ii — t2){coskx — cosky)T^ — 

At 4 sin kx sin kyTi , 



(1) 



where Ti, with i = 1,2,3, are the Pauh matriees and tq is 
the 2x2 identity matrix. The r^ matrices act in orbital 
space. Note that £,xyi}^) must transform as the Aig rep- 
resentation of D4/J; in this representation tq transforms 
as Aig, Ti as 623, and T3 as Big. However, if the degen- 
erate d-orbitals are expressed in terms of the (X, Y) axes 
as [dx, dy), then the tight-binding dispersion becomes)^ 



^xy(k) = [-(ti +t2){cos ka; + cos ky)- 

4^3 COS kx cos ky — ^]tq 

-{h - t2){coskx - cosky)Ti- 

Ati sin kx sin fcj,T3. 



(2) 



Notice that the XY basis is chosen just for the orbitals 
while in real space the system of coordinates is still given 
by {x,y). ^xY(k) also has to transform as Aig which 
means that, when the orbitals are defined in the {X, Y) 
basis, then ti transforms as Big, and T3 as B2g. It 
can be shown that since ri is the matrix that indicates 
inter-orbital electron hopping, this kind of hopping hap- 
pens between nearest-neighbors [ncxt-nearest-neighbors] 
in the (X, Y) [{x, y)] representation. 

As previously discussed, if the values of the param- 
eters are set to <i = —1, ^2 = 1-3, is = ^4 = —0.85 
and fi = 1.54 then the Fermi surface (shown in Fig. [1] 
together with the band dispersion) for the tight-binding 
Hamiltonian is in qualitative agreement with band struc- 
ture calculations for the pnictidcs^ once folding to the 
reduced BZ is performed. Note that the system is half- 
filled (two electrons per Fe site on average) and, due to 
the orbital degeneracy, each orbital is half-filled as well, 
despite the fact that the bands are not equally filled. 
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FIG. 1: (a) Band dispersion and (b) Fermi surface for the half- 
filled two-orbital d-model for the pnictides, with the hopping 
parameters introduced in Ref. |^. 

An important characteristic of the two degenerate d- 
orbitals in this model is that around the hole pockets 
a spinor describing the mixture of orbitals rotates twice 
on encircling these FS's. The inversion and time rever- 
sal symmetry of the twice degenerate d bands ensures 
that at each k point it is possible to choose real spinor 
wavefunctions that arc confined to a plane. The spinor 



has vorticity ±2 around the hole pockets while there is no 
vorticity around the electron pockets<^ As pointed out in 
Ref. |30, this topological characterization of the hole and 
electron pockets is also a characteristic of all the more 
realistic models for the pnictides that include additional 
orbitals. 



B. s- model 

Let us introduce now a two-orbital model with two 
non-degenerate non-hybridized s-like bands, called si(l) 
and 52(2), with dispersion relations given by: 

£,si (k) = 2<i(cos kx + cos ky) + 4^2 cos kx cos ky — fi, (3) 

and 

^S2 (k) = 2^3 (cos kx + cos ky) + At 4 cos kx cos fcj, — /i + A, 

(4) 
where fi is the chemical potential and A is the energy 
difference between the two bands. The dispersions can 
also be written in the basis (si, S2), i-e., (1, 2), using the 
Ti matrices as in the previous case: 



^s(k) = [{ti + <3)(cosfcj, + cos ky)+ 

2(^2 + ^4) cos kx cos ky — fi+ —]to 

+ [{ti - h) {cos kx + cos ky)+ 

2(^2 ~ ^4) COS kx cos ky ]t3. 



(5) 



It is clear that here both tq and T3 transform like 
Aig and for this reason we call this model the s-modcl. 
In Fig. [21 the band dispersion (panel (a)) and the FS 
(panel (b), red circles) are shown for the parameter val- 
ues h = -0.05, i2 = 0.7, ts = -0.1, ti = 0.3, A = 2.8 
and yu = 1.95. The FS of the d-model is also shown (con- 
tinuous black line) for comparison. They are obviously 
very similar, and precisely the goal of this effort is to in- 
vestigate what kind of magnetic and pairing properties 
emerge from these two models that have nearly equal 
Fermi surfaces. 

The hole pockets at the F and M points nest into the 
electron pockets at X and Y, with nesting vectors (0,7r) 
and (tTjO). The system is half-filled but the individual 
bands/orbitals are not. Note that this is the case with 
the orbitals in the multi-orbital systems proposed for the 
pnictides, where nesting occurs between electron and hole 
pockets at the FS but none of the orbitals is exactly half- 
filled4.^ 



C. Coulomb Interaction 

The Coulomb interaction term in both Hamiltonians is 
the usual one, with an on-site intraorbital (interorbital) 




FIG. 2: (color online) (a) Band dispersion and (b) Fermi sur- 
face of the half-filled two-orbital s-model given by Eq. ^ (red 
circles). The continuous line is the FS for the two-orbital d- 
model. 

Coulomb repulsion U [U'), and a Hund coupling J satis- 
fying the relation U'=U — 2J for simplicity, and a pair- 
hopping term with coupling J'=J M' The full interaction 
term is given by 

Hint = U 2^ "i,o,t"-i,o,4. H ^ y ^nj^gni-a 

i,a i.a 

i,a i,a 

(6) 

where dl ^ ^ creates an electron with spin a at site i and 
orbital a = x,y or 1, 2. Si^a ('^i,a) is the spin (electronic 
density) of the orbital a at site i. 



III. MAGNETIC PROPERTIES 

For a single-orbital model, the magnetic structure fac- 
tor is easily defined as 



with 



UJ\ 



^(k)==^e^''---^(r), 



'■) = ]^X1™(')'^(^ + '')' 



(7) 



(8) 



where A^ is the number of sites of the lattice and 

m(i) = rii^i - rii^i = 4,t'^i,t ~ Aidi^i, (9) 

where rn(i) denotes the net magnetization at site i. 

In a multiorbital system the net magnetization at site i 
is obtained in terms of the magnetization of each orbital 
a, and it is given by 

a a 

(10) 



While Eq. [To] characterizes the magnetization that is 
measured in experiments such as neutron scattering, it is 
natural to define generalized magnetic moments mab{i)^ 
given by 



rriabii) = 4a.t^i-''.t - dl^^^dib.i- 



(11) 



With this definition, a generalized form of the magnetic 
correlation functions will depend on 4 orbital indices: 

uJabcd{r) = j^^ ■mab{i)mcd{i + r). (12) 



Thus, it is possible to define orbital dependent magnetic 
structure factors given by: 



SabcdCk) =^e'^-''ujabcd{r). 



(13) 



These orbital-dependent operators may arise from pro- 
cesses as those depicted in panel (a) of Fig. [U where 
having different orbitals at the two vertices is possible if 
the orbitals strongly hybridize to form a bandJ^ 



(a) 



(b) 





FIG. 3: (a) Electronic process that gives rise to the orbital 
components of the structure factor, (b) Same as (a) but in 
the band representation. 

The total orbital magnetic structure factor can then 
be defined as: 



^To(k) 



a.b.c.d 



5*a6cd(k). 



(14) 



Note that there are M'^ orbital dependent components 
of the generalized magnetic structure factor, where M is 
the number of active orbitals in the system. The mag- 
netization that is measured in neutron scattering exper- 
iments is given by Eq. 1101 which in terms of the compo- 
nents of the tensor mab becomes 



m(i) 



/^ 



ni,a,-\~ni^a,i = ^maa(i) = tr[mab{S)]- (15) 



Since m{i) is a trace its value is independent of the basis 
chosen to define the orbitals and it allows to calculate 
the experimentally measured local magnetization. 

Notice that m{\) is the operator that has to be con- 
sidered in order to construct the so-called homogeneous 



or diagonal structure factor defined in terms of the diag- 
onal (intra-orbital) magnetic moments TJiaai'^) and given 



S 



MO 



«4i: 
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a.b.r.i 



- Z^Saa 
a,b 



"laa(i)"l6fc(i + r) = y] SaabbCk.). 

(16) 

Smo is the physical magnetic structure factor that has 
to be calculated in the context of multiorbital systems 
to compare with neutron scattering results.— ^^^ Several 
authors have pointed out the existence of the generalized 
components of the magnetic susceptibility both in the or- 
bital representatio n^ -'^ i ^^ and in the band representation^ 
It has also been pointed out that an orbital-transverse 
density-wave (OTDW) ordered state characterized by 
the non-homogeneous components of the magnetization 
tensor may develop in multi-orbital systems;^ an issue 
that will be further explored and discussed in the present 
work. 



A. Non-interacting case 

In order to understand the relationship between StOi 
Smoi E^nd the properties of the FS of the system, it is 
illuminating to consider the non-interacting case which 
can be easily studied in momentum space. Via a Fourier 
transform of dl^ ^ and di^a,a; Sabcd in Eq.[I3]can be writ- 
ten as 

Sabcdi^) ^ 2_^ (-1)'^^'^ rfq^Q^crC^q+k,6,CTrfp,c,cr''^P-k,<i.cr'- 
p.q.(T,(T' 

(17) 
In momentum space it is natural to use the band repre- 
sentation in which 



PiqiCr.cr',/^.,/^/,^/,!/' 
< ^|a >q< b\^J.' >q+k< 1^\C >p< d\u' >p_k 



(18) 



di 



7t 



'q,Ai,cr'^q+k,A'',cr"p,,y,cr'C'p-k,l/',o-' 



where d'^ ^ ^ creates an electron with momentum p and 
z-spin component a at band i^, while < v\a >p is the 
matrix clement for the transformation from orbital to 
band representation. 

In the band representation, the electronic processes 
that contribute to the magnetic correlations are shown 
in panel (b) of Fig. [3] Since the electronic band cannot 
change as the electron created at the right vertex is de- 
stroyed at the left vertex, in the band representation we 
can define band-dependent components of the structure 
factor given by 

i) ^luv ^i\)^) = / ^ "q,^,CT^q+k,l',o-"p,,y_cr'^p-k,A',o-5 (.lyj 



where the greek indices label the bands. A total structure 



factor can be defined in terms of S^i^i,^ as 

5Ti3(k)=^V-^(k). (20) 



fi,i^ 



Also the homogeneous or diagonal magnetic structure 
factor Smb , analogous of Smo j can be defined as 



Smb (k) — 2_^ '^MMMM 0^) ' 



(21) 



since in the band representation S^^m, — 0, ii ^ ^ v. 
Note that the band representation is the natural starting 
point in approaches based on fermiologyi^'^ 

In the noninteracting case being considered in this sec- 
tion, it is easy to show that 

5^..M (k) = 2 ^ /^ (q) [1 - /, (q + k)] , (22) 



where /^(q) is the Fermi function for the band /i. We 
also find that the components of the structure factor in 
the orbital representation are given by 

S'afccd(k) = 2 ^ < ^i\a >q< b\v >q+k 

q,/J,l> (23) 

< V\C >q+k< d\^i >q /^(q)[l - JM + k)]. 

From the expressions in Eqs. [22] and [23l it can be shown 
that Sto = Stb and Smo ~ Smb only if the orbitals do 
not hybridize to form the bands, i.e., the matrix elements 
are the elements of the identity matrix. In case of a 
nonzero hybridization, then the structure factors in the 
band and orbital representations are different. 



B. d-model 

Numerical Lanczos calculations for the homogeneous 
(or diagonal) magnetic structure factor Smo have al- 
ready been performed in previous literature for the two- 
orbital d-model indicating a tendency towards a magnetic 
stripe ordering for the undoped case, characterized by 
peaks at k = (7i',0) and (0, tt) in SmO'^ This tendency 
is already apparent even in the non-interacting cas o^'^'^^ 
as illustrated in panel (a) of Fig. H] where Smo calculated 
in a 16 X 16 cluster is shown with open circles, along the 
directions (0,0) — (tt, 0) — (7r,7r) — (0,0) in the unfolded 
BZ. The broad peak at k = (7r,0) is clear and it can be 
compared with the curve denoted by the star symbols in 
panel (b) of the same figure where results for the -\/8 x -y/8 
cluster that can be studied numerically exactly (with the 
Lanczos algorithm and for any value of the Hubbard cou- 
plings) are presented. This same behavior is also appar- 
ent in the total orbital structure factor Sto (k) indicated 
by the diamonds in Fig. |4l^a). 

On the other hand, a calculation of the magnetic struc- 
ture factor using the band representation, i.e. S'Ms(k) 
indicated by the squares in panel (a) of Fig. SI shows a 
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FIG. 4: (color online) (a) Magnetic structure factors, total 
and homogeneous as indicated, for the non-interacting two- 
orbital d-model on a 16 x 16 lattice, (b) Homogeneous or- 
bital magnetic structure factor SMo(k.) for the interacting 
case with J/U — 0.25 and at the indicated values of U. The 
results were obtained numerically using an 8-sites cluster and 
the Lanczos method, (c) Total orbital magnetic structure fac- 
tor 5'To(k) for the interacting case, for the same parameters 
and technique as used in (b). 



rather different behavior: instead of a clear peak at (tt, 0) 
there is a featureless plateau around (tt, 0) that extends to 
(7r/2,7r/2). This example demonstrates the importance 
of the matrix elements in Eq. [TH] which differentiate be- 
tween Smo and Smb- In the non-interacting case, both 
functions can be expressed in terms of the Fermi func- 
tions as in Eqs. [22] and [23] allowing us to conclude that 
the peak at (tt, 0) arises from the matrix elements rather 
than from purely nesting effects of the Fermi surfaces. 
Ignoring the matrix elements, it is interesting to note 
that a feature at (tt, 0) can also develop if all the compo- 
nents of the structure factor in the band representation 
are considered and S'TB(k) is calculated, as shown by the 
curve indicated with triangles in Fig. [4l^a). 

The contribution of the band- and orbital-resolved 
components of the structure factor in the non-interacting 
case are presented in panel (a) of Fig. [S] The components 
of the structure factor that contribute to Smo arc Saabb 
with a (b) taking the values x (y) and y (x) indicated 
by the diamonds in the figure, and Saaaa (indicated by 
the circles and squares). It is clear from the figure that 
the peak at (tt, 0) in Smo at the non-interacting level is 
mostly due to the Saabb that arise from the nesting of the 
two bands that contain the same orbital flavors due to 
hybridization, while the components of the form Saaaa 
show features also at (tt, tt) since this wave vector also 
nests the hole (electron) FS's at F and M {X and Y). 
It can be seen that the non-homogeneous components of 
the form Sabab (diamonds) behave as Saabb in the non- 
interacting case and contribute to form the peak at (tt, 0) 
in the total structure factor Sto [triangles in Fig. [31Ja)]. 
For completeness in Fig. [SJa) orbital resolved structure 
factors of the form Sabba (up triangles) and Sabbb (down 
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FIG. 5: (color online) (a) Orbital-resolved components of the 
magnetic structure factor for the non-interacting two-orbital 
d-model on a 16 x 16 lattice, (b) Band/orbital resolved compo- 
nents of the magnetic structure factor for the non-interacting 
two-orbital s-model also on a 16 x 16 lattice. The index 1 (2) 
labels the lower (upper) band. 

triangles) are also shown; Sabba increase the value of Sto 
at (tt, 0) while Sabbb provide a small negative contribution 
to Sto along the diagonal direction of the BZ. Similar re- 
sults were obtained for all the correlations in which three 
of the four indices arc the same. 

In non-interacting single-orbital systems, as studied for 
the cuprates, the spin and charge susceptibilities have the 
same form for all values of non-zero momenta, and any 
features in these functions arise from the nesting proper- 
ties of the Fermi surface. Naively, the same is expected 
in the case of multi-orbital models but, as it will be dis- 
cussed below, the hybridization of the orbitals plays a 
crucial role. In the d-modcl, the peaks in Smo appear to 
be associated with the nesting of the hole- and electron- 
like Fermi surfaces. In the weak coupling picture, it is 
expected that magnetic order with Q equal to the nesting 
moments stabilizes when repulsive Coulomb interactions 
are added. Our Lanczos calculations for Smo and Sto, 
in panels (b) and (c) of Fig. [3] show that this is indeed 
the case. 

The Lanczos calculated orbital magnetic structure fac- 
tor 5'Mo(k), using a -s/S x \/8 sites cluster, is shown in 
Fig. \Mh) for different values of U and at J/U = 0.25. 
This structure factor has a peak at k = (tt, 0) (and 
(0,7r) as well, not shown) that becomes sharper as U in- 
creases, indicating a tendency towards robust magnetic 
order. Mean-field calculations based on these results, but 
extended to much larger systems, indicate that actual 
magnetic order develops at a finite value of UJ^^ 

The Lanczos-evaluatcd behavior of the SmoO^-) peak 
at k = (tt, 0), as a function of U, is shown in Fig.[ni[a), for 
two different values of J/U (0.05 and 0.25). The tendency 
towards a robust magnetic state with increasing U and 
J/U is again clear. 

As previously stated, Smo is the magnetic structure 
factor calculated in the literature for comparison with ex- 




FIG. 6: (color online) Orbital magnetic structure factor at 
wave vector (7r,0) calculated numerically (Lanczos). (a) Re- 
sults for the two-orbital d-model, as a function of the Coulomb 
repulsion U and for the values of J/U indicated, (b) Same as 
(a) but for the s-model. 



periments, but for completeness and for the sake of com- 
parison with the s-model results, in panel (c) of Fig. 21 
we present the Lanczos calculated values for the total 
generalized magnetic moment Sto for the d-model as a 
function of C/, for the case J/U = 0.25. It is clear that 
for the d- model Sto niimics the behavior of Smo- An 
important question to ask is what are the components of 
the orbital-resolved magnetic structure factor that drive 
the development of a peak at Q = (7r,0) (and (0,7r)) 
when the Coulomb interactions are active. In Fig. [7] par- 
tial sums over selected components of the structure factor 
are shown with summations performed over repeated in- 
dices. In panel (a) of Fig. [7] it can be clearly observed 
that Saabbi whose sum over a and b arc indicated by the 
plus signs and the continuous lines in different shades for 
the different values of the interaction, are the components 
that drive that magnetic behavior. In fact, these are the 
homogeneous components that contribute to the physical 
magnetic structure factor Smo- It is interesting to note 
that while J^a.b ^aabb is equal to J2a.b ^abab in the non- 
interacting system (panel (a) of Fig. [5]) the partial sum of 
the non-homogeneous component J2a b ^abab [x symbols 
and dotted lines in Fig.[71[a)] does not increase with U at 
Q while the partial sum ^^ j^ Saabb clearly does. 



C. s-model 

Let us now carry out a similar analysis but for the two- 
orbital s-model defined by Eq.[5] Since in this model each 
band is defined by a single orbital, then it is clear that 
Smo = Smb and Sto = Stb ^ Note that studies based 
on fermiology assume that if hole and electron FS's are 
nested via a momentum vector Q, then spin density wave 
order will arise from a logarithmic instability that devel- 
ops in the spin response at Q and is stabilized by the 
Coulomb intcractioui^i^ In this scenario the spin-density 
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FIG. 7: (color online) Orbital-resolved components of the to- 
tal structure factor (sums over repeated indices are implied): 
Saabb (plus, full line), Sabab (x, dottcd line) and Sabba (star, 
dashed line), for the values of U indicated, obtained numeri- 
cally (Lanczos) at J/U — 0.25 using an 8-sites cluster for (a) 
the d-model and (b) the s-model. 



wave originates from the formation of particle-hole pairs, 
excitons, belonging to the electron and hole FS's (exci- 
tonic mechanism)]^ Our goal is to investigate whether 
this mechanism is valid for the s-model. 
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FIG. 8: (color online) (a) Magnetic structure factors (total 
and homogeneous) as indicated for the non-interacting two- 
orbital s-model using a 16 x 16 lattice, (b) Homogeneous 
orbital/band magnetic structure factor S'A/o(k) for the in- 
teracting case with J/U = 0.25, at the indicated values of 
U. The results were obtained numerically via the Lanczos 
method using an 8-sites cluster, (c) Total orbital/band mag- 
netic structure factor Sto (k) for the interacting case with the 
same parameters as in (b). 

The magnetic structure factor Smo in the non- 
interacting limit, denoted by the squares in panel (a) 
of Fig. m does not show the features expected from the 
nesting of the two Fermi surfaces at momentum Q. The 
structure factor is actually rather flat on all the BZ, van- 
ishing at k = (0, 0) and (tt, tt). These results are not what 



it would have been expected from the nesting properties. 
Note that the results for Smo in the non-interacting 
s-model [squares in Fig. [HIa)] are actually identical to 
the results for the homogeneous structure factor in the 
d-modcl in the band representation Sj^ib [indicated by 
squares in Fig.|4l^a)] , since both systems do have the same 
FS. However, note how different are the results for the 
d-model in the orbital representation [indicated by circles 
in Fig. HI a)]. This is due to the effect of the matrix ele- 
ments that result from the hybridization of the orbitals, 
which play a crucial role in the magnetic properties of 
the system. This effect can be more clearly appreciated 
when the interactions are added. The behavior of the 
peak in Smo(}^) at k = (tt, 0) was calculated with the 
Lanczos method applied to the s-model, by varying U 
and at different values of J /U using an A^ = 8 sites tilted 
cluster. In Fig. ^ih) it can be observed that for values of 
J < O.IC/ the peak in Smo eventually vanishes. On the 
other hand, for J > O.IC/ a rapid increase in the peak's 
magnitude suddenly occurs at a value of U that decreases 
as J/U increases. The increase of the peak at (tt, 0) with 
increasing U is contrasted with the behavior of the fea- 
ture at (7r/2,7r/2) displayed in Fig. |8ljb). Examination 
of the numerical (Lanczos) ground state indicates that 
at this point the Hubbard interaction is strong enough 
to hybridize the two bands and develop magnetic stripe 
order. 
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FIG. 9: (color online) Mean-field calculated orbital/band re- 
solved magnetic order parameters for the s-model, as a func- 
tion of U and for the indicated values of J/U . (a) Total 
homogeneous magnetic order parameter m = m\i -{- 7TI22; (b) 
mn; (c) 77122; (d) mi2 = m2i. 

Based on the numerical results discussed above a 
Hartrec-Fock mean-field calculation was performed, fol- 
lowing technical aspects already widely discussed in pre- 
vious literature.— By this procedure we found that the 
total (homogeneous) magnetization m shown in panel (a) 
of Fig. [S] mimics the behavior of 5Mo(7r, 0). Here, the 
transition to the magnetically ordered state is very rapid, 
resembling a first-order transition. We observed that MF 
magnetic order develops only if J > O.IC/ which is in 
agreement with the Lanczos results shown in Fig. [BJb). 



The mean-field results also indicate that a full gap char- 
acterizes the magnetic state which is then an insulator 
as it can be seen from the MF calculated spectral func- 
tions A(k,uj) displayed on panel (a) of Fig.[TUl It is clear 
that the hybridization of the original bands/orbitals due 
to the Coulomb interaction is very strong and the band 
structure has been totally reconstructed. This behavior 
can be understood in the real-space representation. In 
order to develop magnetic stripes in the half-filled sys- 
tem, it is necessary to have a net magnetic moment on 
each site. In the d-model, each orbital is half filled and 
thus contains a spin- 1/2 that can easily be polarized by 
the interaction. In the s-model, on the other hand, the 
orbitals correspond to the bands, and one orbital is thus 
almost filled while the other is almost empty. Then, there 
are far fewer magnetic moments that can be polarized. 

Thus, we observe that in the s-model the peak at Q in 
the magnetic structure factor does not develop from the 
nesting of the FS but from the Coulomb interaction, and 
it occurs fairly suddenly and at a robust value of C/ > 4 
for the hopping parameters used here. Thus, while nest- 
ing appears to be a needed condition for the development 
of the peak in the magnetic structure factor, it is not a 
sufficient condition. The hybridization of the orbitals 
needs to be present such that the matrix elements allow 
the peak to emerge at sufficiently strong coupling. In 
fact, it is necessary that the bands that are connected 
by the nesting vector Q share the same orbital fiavor. 
If this occurs via hybridization, magnetic order can de- 
velop at relatively weak coupling, but if this is not the 
case, the Coulomb interaction would induce magnetic or- 
der only in the strong coupling regime, as we have verified 
by studying the s-modcl. In this case, the magnetic tran- 
sition is also a metal-insulator transition, as observed at 
least within the mean-field approximation. The d-orbital 
model, on the other hand, is known to display an in- 
termediate metallic magnetic phaseJ^ Thus, the present 
results indicate that the s and d models develop simi- 
lar magnetic behavior only in the strong coupling regime 
while in weak coupling, despite the nearly identical Fermi 
surfaces, both models have quite different ground states. 
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FIG. 10; (color onfine) Intensity of the mean-field calculated 
spectral functions j4(k, tj) as a function oi uj — jj, and k for 
the s-model: (a) in the stripe magnetic phase for U = 5 and 
J/U — 0.25; (b) in the phase with orbital-transverse spin 
order for U = 2.5 and J/t/=0.25. 



1. Orbital-transverse spin order 

While the analysis of the results for the s-model pre- 
sented above indicates that, despite the nesting of the 
electron and hole FS's, no magnetic order, as defined 
by the homogeneous operator, develops in weak cou- 
pling, it is instructive to analyze the behavior of the non- 
homogeneous components and the total magnetic struc- 
ture factor Sto- The non- interacting values of Sto on 
a 16 X 16 lattice are indicated by the triangles in panel 
(a) of Fig. [51 There is a feature at (tt, 0) arising from the 
contribution of the interband components of the form 
Sabba = '5'pi/iy^i, shown by the triangles and diamonds in 
panel (b) of Fig. [5] These are the components of Sabcd 
that contribute to the development of the maximum at 
Q = (tt, 0) (and (0,7r)) because the nesting at Q is be- 
tween FS's defined by different bands. However, this 
type of terms are not part of the definition of the ho- 
mogeneous structure factor Smb- On the other hand, 
the components of the form Saaaa indicated with circles 
and squares in Fig.[SJb), have a very flat shape in all the 
BZ and do not produce a sharp feature. Any other com- 
bination of orbital indices does not contribute to Sto as 
shown in Eq. [221 

The effect of the Coulomb interactions on the feature 
at (tt, 0) in Sto has been obtained with Lanczos calcula- 
tions and it can be seen in panel (c) of Fig. [51 The peak 
slowly increases as U raises from to 4. Notice that for 
the same range of values of U the peak in Smo shown 
in panel (b) of the figure does not change. The obvious 
question is whether this behavior indicates a novel kind 
of order in multi-orbital systems. The answer is provided 
via our MF approach that allows us to evaluate the com- 
ponents of the magnetization iJiab- The homogeneous 
magnetization m displayed in panel (a) of Fig. [51 is ob- 
tained as the sum of the intraorbital magnetizations mn 
and 771,22 shown in panels (b) and (c) of the figure. In- 
terestingly, we found that the non-diagonal components 
mi2 = m.21 develop finite values while the diagonal com- 
ponents are zero for values of J/U > 0.1 as shown in 
panel (d) of the figure. At the MF level we can study the 
real space configuration associated to this finite order 
parameter. We have observed that the orbital spins are 
disordered, which is expected by the lack of features in 
Smo{^), but there are ordered generalized spins Gab(i) 
defined as 



close to each other and the actual ground state depends 
on the parameters 1^ 



Gab(i) = d.^^a,a^a,pdi^b,P, 



(24) 



where S arc the Pauli matrices and the orbital indices 
a ^ h. In Fig. [Tl] we show two configurations of Gi2(i) 
that provide the MF ground state associated with the 
peak in Sto at (tt, 0) (and (0, tt)) when mi2 is finite. 
Panel (a) shows a flux configuration that generates peaks 
at (tt, 0) and (0, tt) in Sto and panel (b) shows a stripe 
configuration that produces a peak at (0, tt). The peak at 
(tt, 0) is generated by a companion configuration rotated 
by 7r/2. Flux and stripe configurations have energies very 
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FIG. 11; Schematic representation of the real space mean-field 
calculated ground states for the s-model when 77112 is non- 
zero, (a) Flux phase; (b) Stripe phase. The dots indicate the 
sites and the arrows represent the MF value of the generalized 
spin Gi2(i) defined in the text. 

The new phase hinted at by the Lanczos calculations 
and stabilized in the MF calculations is insulating. The 
MF calculated spectral functions are shown in panel (b) 
of Fig. [TUl A full gap has developed at the FS indicat- 
ing that this order, if realized, would be observed with 
ARPES measurements. On the other hand, neutron scat- 
tering experiments would not detect it. This can be seen 
by performing a rotation in orbital space given by^ 



d] 



,±,. - ^(4i. 



±42,J- 



(25) 



In this new basis the schematic representations of the 
spins are shown in Fig.[T21 It is clear that while the homo- 
geneous spins in the orbitals -t- (black dots) and - (white 
dots) are ordered, the net spin at each site is and thus, 
neutron scattering experiments will not detect the or- 
der because there is no finite local magnetization. These 
phases appear to be a realization of the orbital-transverse 
density-wave (OTDW) order proposed in Ref. [SJ- 
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FIG. 12: Schematic representation of the real space mean- 
field calculated ground states for the s-model when m\2 is 
non-zero: (a) Flux phase; (b) Magnetic stripe phase. The 
black and white dots represent the orbitals -f and - at each 
site and the continuous and dashed arrows represent the MF 
value of the spin at each orbital. 

Summarizing, a careful analysis of the small-cluster 
ground states obtained via Lanczos techniques, and with 
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mean-field approximations in larger clusters, highlights 
the important role that the orbital composition plays in 
the development of magnetic order. 

For the s-model, it is illuminating to consider the be- 
havior of the total magnetic structure factor Sto, see 
panel (c) of Fig. |S1 calculated numerically as the inter- 
actions are increased. There is a weak increase of Sto 
at Q before the sudden jump at C/ = 4. The behav- 
ior of the resolved components displayed in panel (b) of 
Fig. [7] shows that for < [/ < 4 the partial sum over a 
and b of the non-homogeneous components Satab (x sym- 
bols and dotted line) and Sabba (star symbols and dashed 
lines) increases in value at (tt, 0) indicating the stabiliza- 
tion of the orbital-transverse spin phase. For [/ > 4 a 
sudden increase of the sum of the homogeneous compo- 
nents Saabb (plus symbols and continuous line) develops, 
the non-homogeneous components start to decrease and 
homogeneous magnetic order is established. 



D. Weak Coupling: RPA Analysis 

Additional insight into the weak coupling behavior of 
the d- and s-models can be obtained via the diagram- 
matic RPA method. Using this technique, the magnetic 
susceptibility Xabcd(k, iw) was calculated fiiii^ and the 
static structure factor was obtained by integrating the re- 
sults over w.— In panel (a) of Fig. [131 the RPA-calculated 
diagonal or homogeneous structure factor for the d-model 
is presented. The non-interacting result (in agreement 
with the results indicated by the circles in panel (a) of 
Fig. 2]) are denoted by the dashed line, while results at 
U = 2.64, the coupling strength where divergent behav- 
ior is about to occur for the case J/t/=0.25, are indicated 
by the continuous line. In these results the peak at (tt, 0) 
is very prominent both with and without the Hubbard 
interaction on. 

The same calculation performed for the s-model, pre- 
sented in panel (b) of Fig. [T31 gives rather different 
results. The flat behavior in the noninteracting case 
(dashed line), in agreement with the curve indicated by 
squares in Fig. [Sl^a), is replaced within RPA by a curve 
(continuous line) that develops weak features at incom- 
mensurate values of the momentum. Note that there 
were no precursors of these features in the non- interacting 
limit. Eventually the peak the closest to the F point 
along the diagonal direction of the BZ, indicated with an 
arrow in the figure, was found to diverge when U becomes 
larger than 2.67 for J/U = 0.25. This appears to be an 
illustration of a case in which RPA calculations indicate 
magnetic behavior that is unrelated to nesting proper- 
ties. The RPA results show that an excitonic weak- 
coupling picture in which magnetic order characterized 
by the nesting momentum Q is expected to occur can be 
misleading if the orbital composition of the bands is not 
incorporated into the discussion. In the excitonic picture, 
the expectation is that the Coulomb interaction will al- 
low the formation of electron-hole pairs with the electron 



(hole) in the electron (hole) Fermi surface. Since Smo 
incorporates intraorbital electron-hole pairs, an RPA re- 
sponse requires that the nesting vector connects parts 
of the electron and hole bands that contain the same or- 
bital flavor. This is the case in the d-model where even in 
the weak coupling regime the (tt, 0) magnetic-stripe state 
with two electrons with parallel spins at every site of the 
lattice has the largest weight in the ground state accord- 
ing to our Lanczos numerical studies. Since both orbitals 
are degenerate, the energetic penalization for populating 
both orbitals is U' and there is a gain given by J if both 
spins are parallel. As discussed before, in the s-model, 
on the other hand, the orbitals are non-degenerate and, 
thus, in addition to U' there is an energy A of penaliza- 
tion when two electrons are located in different orbitals 
at the same site. This energy can be larger than the gain 
obtained from J by having parallel spins or than the U 
penalization that arises from introducing both electrons 
in the same orbital. Then, a magnetic "stripe" state can 
only develop when U is comparable to the splitting A. 
This regime, which develops in strong coupling accord- 
ing to our Lanczos and MF calculations, is not captured 
by the weak-coupling RPA method. However, it will 
be shown that RPA is effective at finding the orbital- 
transverse spin state presented in the previous section if 
the generalized structure factor Sto is calculated. 
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FIG. 13: (color online) RPA calculated magnetic structure 
factors for J/U — 0.25, at the values of U indicated (full line). 
The non-interacting results are indicated with dashed lines. 

(a) Homogeneous magnetic structure factor in the d-model. 

(b) Same as (a), but for the s-model. The arrow indicates 
the peak that grows the most as the critical U is reached, (c) 
Generalized magnetic structure factor for the d-model. (d) 
Same as (c), but for the s-model. 

The values of Sto obtained with RPA are presented 
in panels (c) and (d) of Fig. [13] for the d and s models, 
respectively. Both develop a peak at the nesting wavevec- 
tor. The generalized structure factor takes into account 
electron-hole pairs formed by an electron and a hole in 
different bands that are allowed to have different orbital 
flavors. This is the reason why a peak develops now in 
both cases. While in the case of the d-modcl the behavior 
of Sto mimics Si\io and the divergence in both occurs at 
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the same value of U (slightly above 2.64 for J/U = 0.25) 
indicating that the stripe-magnetic order is the cause, in 
the s-model the peak in Sto develops at a lower value of 
U {U = 1.72 for J/U = 0.25) and it results from the or- 
dering revealed by the inhomogeneous components 6*1221 
and S'2112 of the structure factor, i.e., orbital-transverse 
spin order as discussed in the previous section. In this 
new light, we see that the divergence in Smo should be 
disregarded since it occurs for a much larger value of U 
than the divergence in Sto ■ These results show that if all 
the elements of the susceptibility tensor are considered, 
RPA calculations are able to determine the development 
of new ordered phases that can develop in multi-orbital 
systems. Conversely, in multi-orbital systems in which 
orbital-transverse order develop, RPA calculations using 
only the homogeneous susceptibility may lead to unphys- 
ical results. 



E. Strong Coupling Regime 
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FIG. 14: (color online) (a) Fermi surface with its orbital com- 
position for the case of a three-orbital model for the pnictides. 
(b) Mean-field calculated orbital-resolved magnetization for 
the same three-orbital model. The figure was taken from 
Ref. [la for illustration. 



In the regime where the coupling U is sufRciently 
strong such that even in the s-model it is energetically 
favorable to locate two electrons with parallel spins at 
the same site (and in different orbitals), both the s- and 
d-models can be mapped into effective t— J — J' models 
and an insulating state with magnetic stripes can oc- 
cur. In this case the Hubbard repulsion has effectively 
hybridized both bands causing large distortions and ac- 
tually opening a full gap [see Fig. [TUfa)]. In this strong 
coupling regime both models appear to have similar prop- 
erties, but an insulating magnetic behavior does not re- 
produce the experimental behavior observed in several 
of the undoped iron pnictides (such as the 1111 and 122 
families). However, this regime could be applied to the 
chalcogenides: if U is sufficiently strong the magnetic 
behavior that develops in the strong coupling limit is 
more related to the hopping parameters and superex- 
change than to the weak-coupling nesting properties of 
the Fermi surface. While the values of the hopping pa- 
rameters in the Hamiltonian are crucial to achieve nesting 
in weak coupling^S systems in which nesting is not per- 
fect can develop stripe-like magnetic order if they map 
into a t—J—J— model in the strong coupling limit such as 
in the case of the three-orbital model for the pnictides J^ 

The results in this section indicate that in the case of 
the pnictides, even if the five d orbitals are considered, 
the xz and yz orbitals arc the most likely to produce 
the strongest contribution to the metallic stripe magnetic 
order at weak or intermediate values of the Hubbard in- 
teraction because they arc the major constituents of the 
FS's with better nesting and because they are degener- 
ate and, thus, there is no energy A that needs to be 
overcome by the interaction. This is apparent already in 
the three-orbital model for the pnictides, where a mean- 
field calculation shows that magnetic order develops at a 
finite value of U (see panel (b) of Fig.[T3]) J^ In Fig. fTir a) 
it can be observed that the orbital with the best nesting 



associated with Q = (tt, 0) is the yz one, indicated by 
the continuous line. A mean-field calculation of the or- 
bital resolved magnetization maa for a^xz, yz, and xy, 



shows that 



^yz,yz 



grows very rapidly at the lowest value 



of U. The magnetizations for the other orbitals develop 
as U hybridizes and distorts the original bands. Thus, 
in the intermediate U regime when magnetism develops, 
the xz and yz orbitals are the ones that would develop 
the stronger magnetization (albeit for different values of 
Q) giving rise to a magnetic metallic phase. Thus, nest- 
ing seems to drive the magnetization of the xz/yz or- 
bitals while the additional orbital hybridizations that de- 
velop due to the reconstruction of the FS then drives the 
smaller magnetization in the remaining orbitals. 



IV. PAIRING SYMMETRIES 

Regarding the symmetry of the pairing operators cor- 
responding to the models analyzed here, previous nu- 
merical calculations have indicated a competition be- 
tween Aig, B2g, and E^ states in the d-modelfi^i^ as 
shown in panel (a) of Fig. [151 The Eg states corre- 
spond to a p-wave spin-triplet state that becomes desta- 
bilized upon the addition of binding-enhancing Heisen- 
berg tcrmsi^ The favored pairing operators with the 
symmetry Ai^ are all trivial in their orbital composi- 
tion, i.e. they are intra-orbital with the form D^aoD 
where D^ = (dj^. ^ ^, d_j. , ) in the {x, y) basis, and they 
remain intraorbital in the {X, Y) basis. However, the 
B2g pairing operators have a non-trivial orbital composi- 
tion given by D^aiD in the basis {x,y), indicating that 
the pairs arc made of electrons in the two different or- 
bitals. In the {X, Y) basis the B2g pairing operator be- 
comes D'^'asD' with D'l" == (4,^,^ '^-k.r.i)- Thus, in the 
(AT, Y) representation the 62^ pairs are intraorbital but 
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there is an important sign difference between the pairs 
in the different orbitals which makes the orbital contri- 
bution intraorbital but non trivial. It is interesting to 
observe that the intraorbital Big state found with RPA 
calculations in the five-orbital model for the pnictidesii 
would become interorbital in the {X, Y) basis. 




FIG. 15: (color online) (a) Relative symmetry between the 
undoped and the electron-doped ground states for the case of 
the d-orbital model, varying J/U and U. The results were 
obtained numerically via the Lanczos method using a small 
cluster with N = 8 sites (and following steps already dis- 
cussed in previous literature).—— The circles indicate states 
with Eg symmetry, squares correspond to B2g , and diamonds 
represent Aig symmetric states, (b) Same as (a) but for the 
s-model with the triangles denoting Big symmetry. In the re- 
gion above the continuous line the two added electrons form a 
bound state. The dashed line indicates the boundary for the 
stability of the magnetically ordered insulating (MOI) region 
in the undoped state. 

The results for the s-model regarding pairing proper- 
ties are different from those in the d-model. Following 
previous investigations)^ using the Lanczos method we 
have calculated the relative symmetry between the un- 
doped (number of electrons N^ = 16) and electron-doped 
{Ne = 18) states, as an indicator of the possible pairing 
symmetry in the bulk limit. The results are presented in 
panel (b) of Fig. [121 varying U and J/U. For small values 
of U and J/U the doped ground state has symmetry Aig 
in agreement with the d-model, although in a different 
regime of couplings. Increasing U and J/U, the s-model 
ground state switches to the Eg symmetry, i.e. p-wave. 
This p-symmetry arises from the spatial location of the 
electrons since the orbital contribution is trivial. We have 
observed that in the small cluster studied here the spin- 
triplet state with k — (0, 0) is almost degenerate with a 
spin-singlet state with k = (tt, tt). The possibility of hav- 
ing a spin-singlet p-wave state with wavevector k = (tt, tt) 
has been previously discussed long ago in the context of 
the single-orbital Hubbard modeli^ In the present case, 
we need to remember that here k is a pseudo-momentum 
and in the folded representation k = (tt, tt) actually maps 
into (0, 0) so that the actual Cooper pair, if stabilized, 
has zero center-of-mass momentum, but the components 



of the pair belong to bonding and antibonding bands that 
could become hybridized for the large values of the inter- 
actions needed to stabilized these states. As indicated 
in the figure, it was also found that the p-states show 
binding in the small system studied here. In addition a 
small region of bound states with Big symmetry is found 
at even larger couplings. While in the d-modcl our nu- 
merical results indicate that the orbital degree of freedom 
plays a crucial role in the symmetry of the pairing states, 
we observe that this is not the case in the s-model. This 
result seems to indicate that interorbital Cooper pairs are 
likely to be present in multi-orbital systems with strongly 
hybridized bands as it is the case of the pnictides. 

Understanding more deeply why the s-model develops 
its particular pairing properties is at this point unneces- 
sary since the model simply provides an illustration of a 
system with a similar FS as the d-model, and the goal 
of this work was to show that the orbital composition of 
the bands plays a crucial role in determining the symme- 
try of the doped states. The examples that have been 
discussed here clearly show that models with the same 
Fermi surface and the same interactions can have very 
different pairing properties depending on the degree of 
hybridization of the orbitals. It also seems, according 
to the present results, that the relevance of the orbital 
degree of freedom in determining the pairing symmetry 
is also influenced by the degree of hybridization among 
those orbitals. 



V. CONCLUSIONS 

Summarizing, numerical and analytical calculations 
have been performed in order to compare the proper- 
ties of two band models with identical FS's and interac- 
tions, but differing in the degree of hybridization of the 
orbitals to form the bands. Despite the nesting proper- 
ties of the FS's it was discovered that both models have 
similar magnetic (insulating) ground states in the strong 
coupling limit, but they are very different in weak and 
intermediate coupling. The s-modcl offers an example 
in which despite the nesting of the FS and the presence 
of Coulomb interactions, magnetism does not develop in 
weak coupling. However, it was discovered that instead, 
as a result of the nesting in weak coupling, the Coulomb 
interaction stabilizes an orbital-transverse spin ordered 
state with no local magnetization. This state is insulat- 
ing and is characterized by a gap that could be observed 
in ARPES experiments. However, due to the lack of lo- 
cal magnetization, neutron scattering experiments would 
not detect the development of "generalized spin order". 
In fact, standard RPA calculations in the s-model lead to 
incorrect results such as incommensurate magnetic order 
in the physical homogeneous channel. However, when 
the non-homogeneous components of the susceptibility 
arc taken into account, RPA reveals the existence of the 
orbital-transverse spin phase for values of U lower than 
the ones needed to observe the unphysical magnetic state. 
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It is clear that the physical (homogeneous) magnetic 
structure factor depends strongly on the orbital flavor 
of the bands and for this quantity to develop a peak in 
weak coupling it is necessary that the portions of the FS 
connected by nesting have the same orbital flavor. 

The possibility of "hidden" magnetic ordering in the 
pnictidcs has been proposed by several authors^Srd^ as an 
explanation for the unexpectedly low value of the magne- 
tization in several of these materials. The hidden order 
proposed by these authors was "diagonal" , as the config- 
urations we presented in Fig. [T^] after transforming our 
non-diagonal results into a rotated orbital basis. How- 
ever, in multi-orbital systems with more than two or- 
bitals, it may be necessary to consider the non-diagonal 
order as well. In theoretical and analytical calculations 
these non-diagonal hidden orders are revealed by consid- 
ering all the components, homogeneous and inhomoge- 
neous, of the magnetic susceptibility. On the experimen- 
tal side, ARPES can detect gaps that are opened due to 
the "hidden" magnetic order but the traditionally used 
techniques to detect homogeneous magnetic order, such 
as neutron scattering, will fail due to the lack of a local 
magnetization. 

We also found indications of quenching of the orbital 
degree of freedom in systems with non-hybridized or- 
bitals. The orbitals do not appear to play a role in deter- 
mining the symmetry of the pairing states. This degree of 
freedom, though, is crucial in systems with hybridized or- 
bitals. In the case of the pnictidcs in particular, we have 
shown that the ground states with d symmetry found in 



the literature in models for the pnictidcs, such as the Big, 
can be made interorbital by changing the basis in which 
the degenerate xz and yz orbitals are defined. 

The results provided by this work may explain why 
the end member of the series BaFe2(Asi_a;P2;)2 is non- 
superconducting despite displaying the best nesting of 
all the compounds in the seriesi^ If superconductivity 
necessitates magnetic fluctuations they may not be suffi- 
ciently strong in this compound if there is no good match- 
ing between the flavor of the orbitals in the nested bands. 

Finally, our results confirm the perception expressed in 
the analysis of recent photoemission experiments^^ that 
the weak coupling nesting mechanism would not be ap- 
plicable if indeed a hole-pocket band dominated by the 
orbital 3z^ — r^ (with no nesting partner in the electron- 
pocket band) does develop a robust superconducting gap. 
Confirming and then understanding the results of those 
recent photoemission experiments is indeed very impor- 
tant for the clarification of several intriguing issues in the 
challenging physics of the pnictidcs. 
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